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Abstract 

Surface effects on the electronic energy loss of charged particles en- 
tering a metal surface are investigated within linear-response theory, 
in the framework of time-dependent density functional theory. Inter- 
esting phenomena occur in the loss spectra originated by the boundary 
(bregenzung) effect, which is as a consequence of the orthogonality of 
surface and bulk excitation modes. Our calculations indicate that the 
presence of a non-abrupt electron-density profile at the surface severely 
affects the nature of surface excitations, as deduced from comparison 
with simplified models. 
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1 Introduction 

A moving particle approaching a metal surface losses part of its energy as 
a consequence of the creation of electron-hole pairs and both bulk and sur- 
face collective excitations, i.e., plasmons jy, |2[. A theoretical description 
of these electronic excitations is basic to understand the processes involved 
in several spectroscopies, such as x-ray photoelectron spectroscopy (XPS), 
Auger-electron spectroscopy (AES), and reflection-electron energy- loss spec- 
troscopy (REELS) [Q. Equally, the interaction of moving ions with solids has 
represented an active field of basic and applied physics [[§, !J. 
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In this paper, we follow the theoretical framework developed in Ref. 
and focus on the interaction of charged particles moving with constant ve- 
locity along a definite trajectory that is perpendicular to a jellium surface. 
Previous linear-response calculations of the energy loss of charged particles 
entering a metal surface were carried out with the use of simplified models 
|7|. H . Here we present self-consistent calculations of the energy- loss spectra, 
as obtained in the framework of time-dependent density-functional theory 
(TDDFT) ||. Although in the case of charged particles moving inside a 
solid nonlinear effects may be crucial in the interpretation of energy-loss 
measurements [PH[] , surface effects are expected to be well described within 
linear-response theory, unless the velocity of the probe particle is small com- 
pared to the Fermi velocity of the solid [TT[ . 

Contributions to the energy loss coming from the excitation of bulk and 
surface plasmons are studied in detail, as well as the influence of an appropri- 
ate characterization of the electron density at the surface. The latter is found 
to be essential when the particle penetrates the solid, as the surface region is 
particularly perturbed when the particle crosses the surface. Besides, inter- 
ference effects occurring between the various parts of the ion trajectory are 
more relevant than in the case of charged particles moving at a fixed distance 
from the surface, which claims for the need of a self-consistent treatment of 
the surface electronic response. 

In Section 2, we briefly discuss general expressions for the energy- loss 
spectra of charged particles entering a solid surface along a trajectory that 
is perpendicular to the surface. In Section 3, we report the results of our 
full self-consistent calculations, and the effect that the electronic selvage at 
a metal surface has on the energy-loss spectra is discussed by comparing our 
full calculations with those obtained for electron densities that drop abruptly 
to zero at the surface. Unless otherwise is stated, atomic units are used 
throughout, i.e., e 2 = Ti = m e = 1. 



2 Theory 

Let us consider a recoiless particle of charge Z\ moving with non-relativistic 
velocity v along a definite trajectory that is perpendicular to a metal surface 
(see Fig. 1). The moving particle will be characterized by a classical charge 
distribution. The solid will be described by a bounded free-electron gas 
normal to the z axis, consisting of a fixed uniform positive background plus 
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a neutralizing cloud of interacting electrons of density n(z). 

Within linear-response theory, the probability per unit time and unit 
energy for the probe particle at z to transfer energy uj to the solid is found 
to be § 

72 /.+oo p 

P z (uj) = -—— dz cos[u{z - z')/v] / dqqlmW(z, z r ;q,uj), (1) 
2n z v J-00 J 

where q represents a wave vector parallel to the surface and W(z, z'\ q, uj) is 
the so-called screened interaction 



W(z,z';q,uj) = v(z,z';q) + J dz x J dz 2 v(z, z\\ q) x(zi, z 2 ; q,u) v(z 2 , z'\ q), 

(2) 

v(z,z';q) and x{z, z '',q, w ) being two-dimensional Fourier transforms of the 
bare Coulomb potential and the density-response function of the solid |12| 



respectively. The energy that the probe particle at z looses per unit path 
length due to the creation of electronic excitations in the solid is simply [13j 

A VP 1 coo 

-~T = - duuP z {u). (3) 
dz v Jo 

In the framework of TDDFT, the interacting density-response function 
x(z,z';q,u) is fully determined from the eigenfunctions and eigenvalues of 
the Kohn-Sham equation of density-functional theory (DFT) jl4| and the 
exchange-correlation (xc) kernel f xc accounting for short-range xc effects (see, 
e.g., Ref. which in the random-phase approximation (RPA) is taken to 
be zero. To compute x( z , z'; q, uj), we first take a jellium slab of thickness a = 



6 Xf |15| and electron density equal to the average electron density of valence 
electrons in Al (r s = 2.07), and then solve the Kohn-Sham equation of DFT 
in the local-density approximation (LDA) by following the procedure 



described in Ref. |p.7|| . For comparison, we also consider simplified models 
for the screened interaction W(z, z'\ q, uj) of a semi-infinite electron gas, which 
are derived for electron densities that drop abruptly to zero at the surface: 
These are a hydrodynamic model (HDM) |T8|| , and a specular-reflexion model 
(SRM) which expresses the screened interaction in terms of the bulk 
dielectric function e(q,u). 

3 Results 

First of all, we show results for the probablity P z (uj) as a function of z and 
for selected values of the energy u, which we obtain either in the full self- 
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consistent surface RPA, in the HDM, or in the SRM with the bulk RPA 
dielectric function. We set the velocity v — 4i> (i>o = e 2 /h is the Bohr veloc- 
ity) and Z\ = ±1, and our results can then be used for arbitrary values of Zi, 
as the energy-loss probability is within linear-response theory proportional 
to Z\. 

Fig. 2 shows the probability P z (u) for the moving particle to transfer 
either the energy uj = cu p [cu p = ^/Airn, n being the average electron density, 
i.e., 1/n = (4/3)7rrf], corresponding to the excitation of a bulk plasmon, or 
uj = uj s [ uj s = uj p / y/2] , corresponding to the excitation of a surface plasmon. 
When uj = Up [Fig. 2(a)], the energy-loss probability increases as the pro- 
jectile enters the surface and reaches a constant value deep inside the solid. 
At the surface-plasmon energy [Fig. 2(b)], the probability reaches a maxi- 
mum when the probe particle is located near the surface, and diminishes in 
the interior of the solid where only electron-hole pairs can be excited with 

uj = UJ S . 

It is well known that within a classical model consisting of a semi-infinite 
medium of local dielectric function e{uj) the energy-loss probability P z {oj s ) 
would be maximum at z — 0. However, the actual dispersion of the surface 
plasmon, not included in the classical model, shifts the peak position of 
P z (u) s ) towards the interior of the solid. Moreover, the so-called bregenzung 
or boundary effect reduces the coupling to bulk plasmons by the presence 
of surface plasmons. As the electronic selvage is changed from zero (HD 
and SRM) to its actual structure (RPA), the creation of electron-hole pairs 
increases and the bregenzung effect is more pronounced, thereby yielding 
larger surface-plasmon and smaller bulk-plasmon excitation probabilities. 

Fig. 3 exhibits our full surface RPA (solid line), HDM (dashed line), and 
SRM (dashed-dotted line) calculations of the stopping power, as obtained 
from Eq. (3) as a function of z and with v = 4t> . As the projectile is 
entering the surface, our full surface RPA energy-loss calculations show a 
broader and smoother structure, which is mainly due to the presence of a non- 
abrupt electron-density profile at the surface. Before entering the surface, the 
selvage electronic structure yields a larger energy loss than predicted within 
SRM. Inside the solid the actual energy loss is smaller than in the SRM, due 
to the bregenzung effect arising from the orthogonality of the surface and 
bulk plasmon modes. In the interior of the solid, where the electron density 
is constant, both our full surface RPA and SRM calculations coincide with 
the well-known RPA stopping power of a uniform electron gas (horizontal 
dotted line). On the other hand, the HD calculation deep inside the solid 
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coincides with the bulk HD, which differs from the bulk RPA. 

In conclusion, our full self-consistent RPA calculations of the energy loss 
of charged particles entering a metal surface indicate that a fully quantum- 
mechanical treatment of the electronic response of a metal surface is necessary 
to quantitatively describe the effect of the presence of the surface. A system- 
atic self-consistent investigation of both RPA and beyond-RPA energy-loss 
spectra of charged particles moving along arbitrary trajectories is now in 
progress. 

We acknowledge partial support by the Basque Unibertsitate eta Ikerketa 
Saila and the Spanish Ministerio de Education y Cultura. 
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Figure 1: Illustration of an external charged particle impinging perpendicu- 
larly on a metal surface. 



Figure 2: Energy- loss probability P z (u) versus z for a particle of charge 
Z\ = ±1 moving with velocity v — Avq to transfer the energy (a) uo = uo p or 
(b) oj = uj s to the solid, corresponding to the excitation of a bulk and surface 
plasmon, respectively. The electron-density parameter r s is taken to be that 
corresponding to the average density of valence electrons in Al, i.e., r s = 2.07. 
The result of our full self-consistent RPA calculation is represented by a solid 
line. The results obtained with the use of HD and SR models are represented 
by dashed and dashed-dotted lines, respectively. 



Figure 3: Stopping power, as obtained from Eq. (3) as a function of z and 
with v = Av . As in Fig. 2, solid, dashed, and dashed-dotted lines represent 
the result of our full self-consistent RPA calculations and those obtained 
with the use of HD and SR models, respectively. The horizontal dotted line 
represents the result of a bulk RPA calculation. 
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